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Abstract 

A correlated Gaussian random walk(CGRW) model is proposed as a simple model of animal 
dispersal. The general features of CGRW is described. We will discuss how from this single model a 
number of different kinds of correlated random walk can be studied. The CGRWs in one dimension 
is studied in detail and the special limiting cases are discussed where the probability densities are 
found analytically. Numerical simulations are performed and the results are found to be in good 
agreement with the theoretical predictions. Directional persistence in CGRWs is discussed for the one 
dimensional cases. We will show that correlation does not always give rise to directional persistence. 

1 Introduction 

Random walk models have been used to study a large number of phenomena in physics, chemistry and 
biology [U [2] . In biology it has been used to describe various processes occuring at microscopic level 
inside a cell and also at macroscopic level to understand the motion of the animals in its habitat. The 
simple random walk model is useful in describing the purely random motion which does not show any 
correlation. The simplest example would be the diffusion of particles in a media. However, in a large 
number of cases it has been observed that the random motions in nature are correlated [23] [24] ?]. 
The dispersal of animals, birds, insects etc show correlation and persistence. In the correlated random 
walk(CRW) models the directions of the consecutive steps are correlated. CRWs have been studied 
extensively since the work of Taylor [3 , Goldstein [4 and Patlak [5 j. CRWs on lattices have been studied 
by a number of authors (see Ref. [6] [8] E9 HDI HQ EJ). It' has been found that the CRWs are 
direct ionally persistent [5j EI [15] [16]. The motion of insects like beet els, ants, butterflies etc has been 
studied theoretically using CRWs by a number of authors [UJ [T7J [21] . Two dimensional models of CRW 
has been used to describe the motion of amoeba [19] and insects [13] . In these models [19] [13] consecutive 
steps are correlated via the relative angle with step size being a random variable. A similar model with 
constant step size had been studied Concepcion and Argyrakis [22]. Bovet and Benhamou [20] has applied 
the CRW to explain the spatial patterns of search paths of foraging animals. An interesting work on 
the postural sway of human body had been modeled Collins and De Luca [18] using CRW. Application 
of CRW models to animal search strategies has been studied Bartumes et. al. [23]. They propose that 
some animals may have evolved the ability of performing Levy walks as adaptive strategies to face search 
uncertainities. 

Our main aim in this work is to study CRWs which can be applied to the problem of animal dispersal 
in general. In this article we will study the correlated Gaussian random walk (CGRW) model. In Sec. 

2 we shall define the CGRW model in and discuss the general features of the model. We will see how 
the behavior of the random walk depends on the choice of 'strategy' and the 'memory' of the random 
walker. It is found that 'strategy' and 'memory' both play a very crucial role in animal dispersal. We 
will also identify how various parameters of our model are related to the 'strategy' and 'memory' of the 
walker. Further more, we will show that correlation between consecutive step does not necessarily lead 
to directional persistence. In Sec. 3 we discuss the symmetric Gaussian random walk(SGRW) in one 
dimension. The probability densities for the limiting cases are calculated and are compared with the 



* Email address: bagarti@iopb.res.in 



1 



numerical simulations in Sec. 4. In Sec. 5 the directional persistence in these random walks has been 
discussed. Finally we conclude in Sec. 6 with possible application to correlated motion of animals. 



2 The CGRW in R 



A Gaussian random walk in d dimension of length n is defined by 



n 



k=l 



(1) 



where S n G R d is the displacement of the random walker from the origin after n steps and G R d is the 
fcth step taken by the random walker. If the collection of random vectors xi,X2, . . . x n are independent 
and identically distributed Gaussian random variables then Eq. (1) defines a simple Gaussian random 
walk(SGRW) in R d . A CRW can be defined in the following manner. The probability for each step x& 
that the random walker takes in the fcth step depends on a function of the previous step x&_i denoted 
by /x(xfc_i) G R d for all fc = 2,3, ... , n. Note that /i(xq) = fio is itself a random vector with a given 
probability density p(fio) which can be specified by considering the symmetry of the problem. We define 
a CGRW by Eq. (1) and the Gaussian probability density for the step x& 



for all fc = 2, 3, . . . , n. We observe that a more general CRW can be defined with \i depending not only 
on the previous step but all the previous steps. The probability density P(S n ,/x(xo)) is given by 



Note that /i(xq) is the initial condition. In general it is not possible to integrate the right hand side of 
Eq.(3) however we will see that in one dimension for certain choice of the function /i(x) we will obtain 
the probability density P(S n ,/x(xo)) within some approximations. The probability density P(S n ) = 
J P(S n , /i(xo))p(/io)d/io can then be determined from Eq.(3). 

The Eqs. (1) and (2) with a given p(/io) completely define our CGRW in R d . Let us now discuss how 
correlation arises in this random walk. In Fig. (2) we show a schematic of the random walk. The steps 
Xfc-i and Xfc are denoted by solid arrows. Also the function /x(xfc_i) is a vector which is denoted by the 
arrow with dashed line. The probability that the random walker takes a random step x& depends on the 
length |xfc — /x(xfc_i)| and is given by P(xfc)dxfc. The expression in Eq. (2) suggests that for a fixed a 
and a given /i(-), most of the time the random walker will take step that does not deviate very drastically 
from its previous step. This gives rise to a correlation in the random walk. Now suppose that we have 
a sufficiently large such that |/i(x)| <C a for all x G R d then the ratio |x n — /i(x n _i)|/2<7 2 ~ |x n |/2<7 2 so 
that CGRW becomes a SGRW having no correlations. In the other hand if a is small then the CGRW 
will be very highly correlated. We can associate important physical meanings to a and \i. Since the fcth 
random step is affected by the (fc — l)th step through the function \i it can be called the strategy and the 
quantity 1/a acts as a measure of how closely it follows its strategy //. The quantity 1/a is the memory 
the random walker possess. 

Before we go to the next section let us have a quick look at the the kind of CRWs we shall be studying. 
In Fig. 1 we have plotted a collection of typical trajectories of CGRWs in R that start from the origin. 
The length of each random walk is n = 500 and the Fig. 1(a) and (b) show 100 such trajectories. The 
trajectories in Fig. 1(a) shows persistence where as (b) does not. We will discuss how the CGRW in Fig. 
1(b) has almost 'zero directional persistence'. The random walkers here follow two different strategies 
so it is not surprising that a change of strategy can change the nature of the walks. 



A symmetric CGRW in one dimension is defined by Eq. (1) and (2) with d = 1 and p(/io) uniform. We 
call it a symmetric CGRW because all possible /iq can occur with the same probability. Here a uniform 
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3 Symmetric CGRW in R 
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Figure 1: Schematic of CGRW in two dimensions. 
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Figure 2: Typical trajectories of GRWs of length n = 500. (a) Trajectories of CGRWs (in green) with 
correlation depending only on step direction, the trajectories in red are those of SGRWs, the values 
a = 0.5 has been used, (b) Trajectorie of CGRW for a = 2.5, here correlation depends on both direction 
and the size of the steps. 



p(/jl) implies that the walker will take the first step to the left or to the right with equal probability since 
sgn(/io) = ±1 with probability 1/2. From Eq. (3) the probability density P(S n , /i(#o)) for a given fi(xo) 
can be written as 

/elk 
—e- ikS "F n (k,v(xo)), (4) 

where N = (\^27rcr)~ n and the function F n (k, fi(xo)) is 



J exp {^2 ~ ~ — 2(7* 1 ~ + ikx ^j n dx 3i 



F n (k,{i(x )) = J exp |^2^ ^2 Vikxi ) | | dx :) , (5) 

From Eq. (5) we obtain the following recursion relation 

F n (k,fJi(x)) = J exp ^ -\y-^(x)\ 2 +iky ^j Fn _ l{k ^( y )) dy ^ ( 6 ) 

where n > 1 and Fo(k^jj J ) — 1. We shall use Eqs. (4-6) to obtain the probability density P(S n ) for the 
cGRWs that we discuss in the following. 

The choice of the function n(x) decides the properties of the random walks. The simplest case is 
when n(x) = c is a constant. If the c is zero then we have a SGRW which has no correlations between 
consecutive steps. Also for nonzero c we shall again have a SGRW which is baised in one direction 
depending on the sign of the constant c. An interesting long tailed behavior of the probability density 
can arise if fi(x) is diverging at some values of x. A choice of ji(x) = 1/x can give rise to such long tailed 
behavior. There could be a number of CGRWs which may be interesting but most of them would be 
too difficult to understand analytically. Here we discuss two different random walks for which we can 
calculate approximate analytical results that gives us insight into the nature of these correlated random 
walks. First we will discuss a CGRW in which every step that the random walker takes depends only on 
the direction of the previous step. The second type of CGRW that we discuss will not only depends on 
the direction of the previous step but also on the size of the step (see Fig.l). 
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3.1 CGRW with correlation depending on direction of the steps 

We define fi(xk) = Xk where Xk = sgn(xk) and fio — ±1 with probability 1/2. With this setting Eq. 
(1) and (2) in one dimension defines a symmetric CGRW. Substituting = x in Eq. (6) we obtain a 
recursion relation for F n (k,x) which can be written in the following matrix form 

F„ = ^ae-^/'MF^, (7) 

where F n = (F n (&, —1), F n (k, 1)) T and M is a 2 x 2 matrix given by 



^(l-erff " 1 ^ 2 )) e-^fl + erff - 1 ^ 2 )) , 

M = 2 (8) 

1 e ik (l - erf(±±f^)) e ik (l + erf(±±§^)) 1 ' 



Once F n (k,x) is obtained from Eq. (7) we can find the density P(S n ,xo) from Eq. (4). Note that if 
x,y £ R then Re(erf(x + iy)) = — Re(erf(— x + %)) and Im(erf(x + %)) = Im(erf(— x + %)). Setting 
x = l/y/2a, y = ka j \[2 and defining A = Re(erf(x + iy)) and B = Im(erf(x + iy)), the eigenvalues of M 
are 

ei, 2 = (1 + A) cos /c - £ sin fc ± [-4 A + ((1 + A) cos k - B sin k) 2 } * . (9) 
From the recursion relation Eq. (7) we have 

/ 77- \ n/2 i 2 2 

F„ = (-J a n e~ nk a /2 M n F . (10) 
The matrix M™ can be written in terms of the eigenvalues ei^ and M as 

V ei — e 2 / V e 2 — 6 i / 

where I is the 2x2 identity matrix. Substituting Eq.(ll) in Eq.(10) and using Eq. (4) we obtain 

Q„ = 2"" J ^. e -nk^/2-ikS n (12) 

F , 



X 



M-6 2 I \ + ^ / M-eil 

ei - e 2 / 2 I 6 2 - ei 



where Q n = (P(S n , —1), P(S n , 1)) T and since i*b(fc, ±1) = 1 with probability 1/2 the probability density 
^cGRw(#n) = (P(S n ,-l) + P(S n ,l))/2 can be obtained from Eq. (12). Although it is difficult to 
obtain a closed form expression for the integral in the right hand side of Eq. (12), it can always be 
integrated numeically for any values of a and can be compared with the probability density obtained 
from numerical simulations. In Sec. 4 we will use the expression in Eq. (12) to obtain probability 
densities for the limiting cases a«l and a ^> 1. 

3.2 CGRW with correlation depending on both direction and magnitude of 
the steps 

Now we present the second type of random walk in which ji(x) depends also on the magnitude of x. Let 
us define /n(x) by the function 

f x \ = / sgn(x)p, if \x\ < x , , 13) 
^ ) y 0, otherwise, ^ ' 

where p > 0. The CGRW for the above choice of fi(x) has completely different properties from the one 
we discussed in Sec. 3.1. Substituting the above ji(x) in Eq. (6) we obtain the same matrix recursion 
relation Eq. (7) with F n = (F n (fc, — p), F n (k, 0), F n (k,p)) T and M is now a 3 x 3 matrices 

e-^Vi(-p) e- ipk *l) 2 (-p) e~ ipk ^(-p) 
M= ( ^i(O) ^ 2 (0) ^ 3 (0) | (14) 

e^ViO) e ipk i/j 2 (p) e ipk i/j 3 (p) 
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where = erf(^±^^) - erf(^), ^ 2 (p) = 2 - erf( » + % + *° ) + erf( ^-° ) and V> 3 (m) = 

erf ( /i+^T ^ _ er £^M±l^_z^o An expression of the form Eq.(12) can be obtained for this case by di- 
agonalizing M but will be too cumbersome to write it here. Such an expression can however always 
be numerically integrated to compare with the results obtained from numerical sumulations. We will 
use the recursion relation Eq.(7) with the matrix in Eq. (14) while discussing the limiting cases in the 
following section. We will find an approximate solution to the recursion relations and we will see that 
the results agree quite well with the numerical simulations. 



4 Approximate solutions and numerical results 

It is in general difficult to find a closed form expression for P(S n ) using Eq.(12). Here we discuss the 
approximate solutions when a«l and a ^> 1. When a>l we can approximate |x& — /i(x/ e _i)| 2 /2<j 2 ~ 
|x/ c | 2 /2<j 2 and the CGRW approaches SGRW. The probability density is PsGRw(^n) = exp(— 5 2 /2na 2 )/\/2n7r<7 
. On the other hand when P(S n ) for the CGRW shows a bimodal behavior which can be written 

as the sum of two Gaussians centered around n and —n. Note that in the limit k — >> oo and a > 0, 
erf((±l + ika 2 )/V2a)e~ k a I 2 = and for k -> it is erf(±l/\/2a). The significant contribution to the 
integral in Eq. (12) will come from a small neighborhood of k = 0. Expanding erf((±l + ika)/V2a) in 
Taylor series near the origin we have 

erf((±l + ika)/V2a) = erf(±l/>/2<7) + X(k 2 + i2k) + 0(/c 3 ), (15) 

where A = ae~ 1 l 2(7 ' 1 j \f 7 h\ . Further for A <C cr <C 1 we take the limit A — >> and we can write from Eq. 
(15), A = ±1 and B = 0. Substituting these values in Eq. (9) we have ei = 2e~ lk and e 2 = 2e lk . Now 
using Eq.(12) the probability density can be written as 



p 

r approx 



(S n ) = J g e -^n-»fc 2 CT 2 /2 cQsnk (16) 



The probability density P(S n ) for the limiting cases is given by 



v =^exp{-Sl/2na 2 ), if<r>l, 
P, PP .(S n ) = { -^[exp(^#) + exp(^i#)], i 



2 y/27rna 



if (7<1. 



In Fig. (3) we compare the limiting solutions Eq.(17) and the exact results obtained from Eq. (12) with 
the histograms obtained from numerical simulations. 

We turn next to the CGRW described in Sec. 3.2. In the limiting case xq/ct <gil, p/a ^> 1, the matrix 

M in Eq. (14) assumes a simple form with ^ 2 (0) ~ 2[1 - erf(^)], ^i(O) = ^s(O) ~ ^f^r exp(/c 2 a 2 /2), 
-02 (ip) — 2 and ^i(±p) = ^(zbp) ~ 0. With these approximations we can write from Eq. (7), 



F n (k,±p) = v^(7e ±<fe P- fe2ff2 / 2 F n _i(fe,0), (18) 
F n (fc, 0) = vWe- feV2 / 2 (l - erf(-^-))F n _i(fe, 0) + x^F^k, + ^n-iCMl- (19) 

It is straight forward to solve these recurrence relations but one obtains cumbersome expressions since 
the recursion relation for F n (k,0) is second order. In our analysis we will find an approximate solution 
to the Eq. (18-19). In the limit xo — >■ 0, Eq. (19) becomes independent of F n _i(fc, ±p) and Eq. (18) 
becomes F n (k, ±p) = e ±zkp F n (k, 0). An approximate solution of Eqs.(18) and (19) can be obtained if we 
assume F n -i(k,±p) ~ e ±lkp F n -i(k, 0) for xq ^ 0. We have 



F n (k,0) = 



2^a(l -erf(^))e- feV2 / 2 + 2x cosfcp 
v 2<j 



(20) 



The probability density P a pp.(^n) can be written as 



-^app. (*^n) 



67T 



exp(-zfcS n ) ae" u/2 + ^ cos fcp (1 + 2 cos fep). 



(21) 
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Figure 3: A comparison of histograms and probability densities for various cases, (a) Probability densities 
for a = 0.25 shows a bimodality. (b) At a = 10, CGRW SCGRW. (c) For a = 1 with the red curve 
showing PsGRw(Sn)- (d) A comparison of numerically obtained probability density with the exact 
expression Pcgrw(£) from Eq.(12) and the approximate probability density P a pp.(^n)- 



where a = (1 — eT K^§^)) an d ft — • In Eq.(21) a factor 1/3 is multiplyed since Prob(/io £ 

{— p, 0,p}) = 1/3. As in the previous case for large a, cGRW approaches the sGRW we see the same 
behavior here as well, hence P a pp.(^n) = ^SGRw(Sn)- We are rather interested in the other limit which 
shows interesting oscillating features in the probability density (e.g. see Fig. 4 (a)). Integrating Eq. (21) 
gives 

( T \ -(S n -(r-2t)p) 2 ^±1 ( r I l\ -( Sn -(r+l-2t')p) 2 

In Fig. 4 (a) and (b) we can see Eq. (22) has excelent agreement with the numerical results. In the 
next section we shall discuss the directional presistence of the two CGRWs we have discussed above. 
The numerical simulations of the random walks were performed for n = 500 for values of 0.25 < a < 10. 
In Figure 3(c), 3(d) Pcgrw(Su) was calculated numerically using the trapezoidal rule from the integral 
in Eq. (12). Clearly for large a Fig. 3 (b) and 4 (d) shows that the correlations are lost and the 
cGRW approaches sGRW. The approximate results obtained in Eq.(17) for a <C 1 and Eq.(22) agree 
very well with our numerical simulations. Fig. 3 (a) shows the bimodal probability density P app .(Sn) 
at a = 0.25 and in Fig. 4 (a) P a pp.(^n) shows oscillatory pattern where a = 0.5, xo = 0.001 and 
p = 100. The bimodality is a direct consequence of the persistence in random walk. The root mean 
square displacement is ~ n which suggests that the motion is super diffusive. In Fig. 4 (a) the oscillatory 
nature of the probability density tells us that the random walk in the second case consists of clusters 
of walks between long jumps. This type of walk is actually observed in nature. Imagine for example 
certain insects, say butterflies. The random motion of a butterfly which when flying over a tree looking 
for nectar in flowers will consist of hovering over bunch of flowers and a less frequent long flights of a 
typical. This typical length in our model is p and peaks in Fig. 4(a) denotes the clusters of walks. The 
CGRW model in Sec. 3.2 describe the random motion of such a butterfly with only exception that the 
butterfly is myopic. We note that in the real case the butter fy is not myopic and therefore the strategy 
la has to be designed accordingly. 



app 



-Y 

r=0 



2 r ^/2n(n - rja 
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Figure 4: A comparison of histograms and probability density P a pp.(^n) from Eq.(22) with xq = 0.001. 
(a) Probability densities for a = 0.5 shows an oscillatory behavior, (b) at a = 2.5 the probability density 
has no oscillation and (c) for a = 10.0, CGRW ->> SGRW, here P a pp.(Sn) = Psgrw(S u ) (in red), (d) 
hi(S 2 ) (in red) show the approximate theoretical estimate obtained from Eq. (22) with xo = 0.001 and 
a = 0.5, compared with the numerical results, (S 2 ) for SGRW is shown in (blue). 

5 Persistence in CGRW 

Persistence in the motion is an important property of CRWs. Directional persistence is the property, by 
virtue of which, an object in motion prefers to move in the same direction. For the case of a SGRW, 
since steps are uncorrelated it is very unlikely that the random walker will continue to move in the 
same direction for a long time. The probability that the random walker will continue to move in the 
same direction for n consecutive steps is 2~ n . For the CGRW defined in Sec. 3.1 this probability is 
p n = [^(1 + erf(^=^))] n . Although these probabilities are monotonically decreasing with the number of 

steps n for both SGRW and CGRW, we note that p n > 2~ n for all n > and a finite. This suggests 
that the rate of decrease of p n is slower than that of 2~ n . It is this property of the CRWs which is called 
directional persistence. Persistence in random walks can result in amomalous diffusion. From Eq. (17) 
the mean square displacement (S 2 ) — na 2 for a ^> 1 and (S 2 ) — n(n + a 2 ) for a <C 1 which shows 
anomalous behavior for a <C 1. Interestingly in the limit a — >• oo the probability p n — >> 2~ n (see Fig. 5) 
the correlations are lost and we say the SGRW has 'zero persistence'. A random walk with p n > 2~ n is 
a persistent random walk. Recall that we have associated the meaning 'memory' 1 to the quantity 1/cr, 
hence a vanishing L memory 1 leads to a complete loss of persistence. 

The CGRW defined in Sec. 3.2 has a different story to say. From Eq. (22) the mean square 
displacement is found to be (S 2 ) = (a + f3) n ~ 1 [p 2 (2a + (2 + 3n)(3) + 3naa 2 ]/3. This is an approximate 
estimate valid only when xo/cr <C 1 and p/cr ^> 1 (see Fig. 4(d)). For a vanishingly small xo but with a 
finite memory (i. e. a < oo) we have (S 2 ) —> (2p 2 + 3na 2 )/3. For any xq which can be set as close to 
zero as we please, we will always have steps of length ~ p in the /cth step when ever the (k — l)th step has 
length less than x$. Now consider any two such consecutive steps of length ~ p at the /cth and k'th. time 
step where k r > k + 1 > 0. The steps will have length \xj\ > xo for all j = k + 1, . . . , k r — 1. This portion 
of the CGRW is an unbaised random walk as fi(xj) = for all j = k + 1, . . . , k' — 1. We also observe that 
fi(xj) = ±p with probability 1/2 for j = k and k' as the (k — l)th and (k f — l)th steps can be positive or 
negative with equal probability. The probability p n that this CGRW will move in the same direction for 
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n steps will be 2~( n_2 ) x 1/2 2 = 2~ n . The CGRW described in Sec. 3.2 therefore has zero persistence. 
This fact is also true for any finite xo with p/a ^> 1. This suggests that correlation between consecutive 
steps alone does not always lead to directional persistence in random walks. Further, the calculation of 
the probability p n for the real cases of animal dispersal should give a hint of how to design the function 
fi and the constant a. 




Figure 5: The probability p n for values a = 0.25, . . . 10.0, the fast decaying curve implies least persistence. 

Acknowledgements 

I wish to thank Jim Chacko for the useful discussions and suggestions. 

6 Conclusion 

We have studied a model of CGRW in which correlation between steps has been introduced by a prob- 
ability density function of the individual step. Some limiting cases of CGRW in one dimension were 
studied in greater detail where we were able to find the probability densities of displacement S n . We 
note that the models we have studied in Sec. 3.1 and 3.2 are two particular cases of CGRW. It is 
only through experiments that one can get a clue of how to design an appropriate 'strategy' function 
ji. In the model in Sec. 3.1 we studied the CGRW in one dimension in which the consecutive steps 
were correlated only in their directions. We found that this walk is persistent. Next we investigated 
a model in which the correlation depends also on the length of steps. Here we found that the random 
walk trajectories consists of clusters of walk and long jumps. This showed an oscillating behavior in the 
probability density. It was found that this walk did not show directional persistence. It would also be 
interesting to see the effect of time dependent memory in these walks. Further investigations can be 
carried out to see whether new interesting features arise in higher dimensions. We were able to identify 
parameters which could be associated with the 'strategy' and 'memory' of the random walker. In one of 
the cases we saw that correlation brings directional persistence in random walk but it was also shown 
that directional persistence may not be seen for any arbitrary 'strategy' used by the random walker with 
a finite 'memory'. This model has a potential application to a wide range of animal dispersal. 
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